Noname manuscript No. 

(will be inserted by the editor) 



A modification of Einstein-Schrodinger theory that 
contains Einstein-Maxwell- Yang-Mills theory 

J. A. Shifflett 



Nov 9, 2008 

Abstract The Lambda-renormalized Einstein-Schrodinger theory is a mod- 
ification of the original Einstein-Schrodinger theory in which a cosmological 
constant term is added to the Lagrangian, and it has been shown to closely 
approximate Einstein-Maxwell theory. Here we generalize this theory to non- 
Abelian fields by letting the fields be composed of d x d Hermitian matri- 
ces. The resulting theory incorporates the U(l) and SU(d) gauge terms of 
Einstein-Maxwell-Yang-Mills theory, and is invariant under U(l) and SU(d) 
gauge transformations. The special case where symmetric fields are multiples 
of the identity matrix closely approximates Einstein-Maxwell- Yang-Mills the- 
ory in that the extra terms in the field equations are < 10 -13 of the usual 
terms for worst-case fields accessible to measurement. The theory contains a 
symmetric metric and Hermitian vector potential, and is easily coupled to the 
additional fields of Weinberg-Salam theory or flipped SU(5) GUT theory. We 
also consider the case where symmetric fields have small traceless parts, and 
show how this suggests a possible dark matter candidate. 

1 Introduction 

The Einstein-Schrodinger theory is a generalization of vacuum general rela- 
tivity which allows non-symmetric fields. The theory without a cosmological 
constant was first proposed by Einstein and Straus 1,2,3,4,5 . Schrodinger 
later showed that the theory could be derived from a very simple Lagrangian 
density HEUE] if a cosmological constant was included. Einstein and Schrodinger 
suspected that the theory might include electrodynamics, but no Lorentz force 
was found [9 ,10 when using the Einstein-Infeld-Hoffmann (EIH) methodpTj 

E2]. 
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In a previous paper [13] we presented a simple modification of the Einstein- 
Schrodinger theory that contains Einstein-Maxwell theory. The Lorentz force 
definitely results from the EIH method, and in fact the ordinary Lorentz force 
equation results when sources are included. The field equations match the 
ordinary Einstein and Maxwell equations except for extra terms which are 
< 10~ 16 of the usual terms for worst-case field strengths and rates-of-change 
accessible to measurement. An exact electric monopole solution matches the 
Reissner-Nordstrom solution except for additional terms which are < 10~ 65 
of the usual terms for worst-case radii accessible to measurement. An exact 
electromagnetic plane- wave solution is identical to its counterpart in Einstein- 
Maxwell theory. The modification of the original Einstein-Schrodinger theory 
is the addition of a second cosmological term A z g^ u , where g^ is the symmet- 
ric metric. We assume this term is nearly canceled by Schrodinger's "bare" 
cosmological term A^N^, where N^ v is the nonsymmetric fundamental ten- 
sor. The total "physical" cosmological constant A = Ab + A z can then match 
measurements of the accelerating universe. A possible origin of our A z is from 
zero-point fluctuations p!4 l[T5 l lT6l fT7] and Higgs field vacuum energy, although 
we just take A z as given, without regard to its origin. The theory in [13] is 
related to one in [18], but it is roughly the electromagnetic dual of that the- 
ory, it uses a different nonsymmetric Ricci tensor with special transformation 
properties, it allows coupling to additional fields (sources), and it allows A ^ 0. 

Here we generalize the theory in [13] to non-Abelian fields by letting the 
fields be composed of d x d Hermitian matrices [19] . This is done much as it is 
done in [20l[2"T] with Bonnor's theory [22]. The resulting theory incorporates the 
U(l) and SU(d) gauge terms of the Einstein-Maxwell- Yang-Mills Lagrangian, 
and if we assume that symmetric fields are multiples of the identity matrix, we 
get a close approximation to Einstein-Maxwell- Yang-Mills theory. The theory 
can be coupled to additional fields using a symmetric metric <? Mt , and Hermi- 
tian vector potential A^. If we let d = 2 and couple the theory to the Standard 
Model, the U(l) and SU{2) gauge terms are incorporated together with the 
geometry, and the combined theory is invariant under U(l) <£> SU(2) ® SU(3). 
Likewise, if we let d = 5 and couple the theory to flipped SU(5) GUT theory, 
the U(l) and SU(5) gauge terms are incorporated together with the geometry, 
and the theory is invariant under U(i) <g) SU(5). Note that flipped SU(5) GUT 
theory 24,25] avoids the short proton lifetime and other problems of the orig- 
inal SU(5) GUT theory. Assuming that we use the usual fermion and Higgs 
field Lagrangian we will get the usual energy- momentum tensor in the Einstein 
equations, the usual current in Ampere's law, the usual equations of motion 
for fermion and Higgs fields, and the usual mixing and mass acquisition in A^. 

This paper is organized as follows. In Sj2]we discuss the Lagrangian density. 
In $3] we show that the Lagrangian density is real and invariant under U(l) 
and SU(d) gauge transformations. In $4] we consider the special case where 
the symmetric fields are multiples of the identity matrix, and we quantify how 
closely this theory approximates Einstein-Maxwell- Yang-Mills theory. In $5]we 
consider the case where the symmetric fields have small traceless components, 
and show how this suggests a possible dark matter candidate. 



3 



2 The Lagrangian density 

Einstein-Maxwell-Yang-Mills theory can be derived from a Palatini Lagrangian, 
£(r* T ,g pT , A v ) = -j^V^dWRv^r) + [n-2)A h ] 



1 ^ tr{F pa g^gP-F^) 

"^:v-3 — T^d +£m(9i*',-A 1 „ip,<l>- ■ •), (i) 



containing a metric g vjl , connection rj* , and Maxwell- Yang-Mills field tensor 

F Vfl = 2A M + i[A v ,A p ]g c /hV2d. (2) 

Here At, is a cosmological constant and g c is the coupling constant. The vector 
potential A a is composed of dxd Hermitian matrices and can be decomposed 
into a real U(l) gauge vector A° a and d 2 ~l real SU(d) gauge vectors A", 

A v = IA° v +r a A a v . (3) 

Here I is the identity matrix and the generators r Q are dxd matrices with 

[T a ,T b ] = iV2df abc r c , t* = rj, tr(r a ) = 0, tr(r a T b ) = d<$£, (4) 

where the f a bc are totally antisymmetric structure constants. For example, 
with d=2 the r a are the Pauli matrices, f a bc = £abc, and g c — e/ sind w where 6 W 
is the weak mixing angle. The £ m term couples g^u 

and Afj, to additional fields, 
as in Weinberg-Salam theory or flipped SU(5) GUT theory. The symbols ( ) 
and [ ] around indices indicate symmetrization and antisymmetrization, and 
[A,B]=AB-BA. Note that the term tr{F pa g ap g pv F vp ) in (P) contains both 
U(l) and SU(d) gauge terms. Our conventions differ a bit from the usual ones. 
The factors of 2c? in (|i!2p result because our r a are normalized like the identity 
matrix as in ([4} instead of with tr(r a T b ) = S b l /2, so that A and A a are on an 
equal footing in Q. As in [13], a factor of 1/Air in Q]) results because we are 
not using the Heaviside-Lorentz convention. We are also using geometrized 
units with c = G=l instead of natural units with c=h = l. 

The original Einstein- Schrodinger theory allows a nonsymmetric N pv and 
rh. in place of the symmetric g pv and J 1 *, and excludes the tr{F pa g ap, g pv F Vi) ) 
term. Our theory introduces an additional cosmological term %A Z as in [13j . 
and also allows F? and N vp _ to have dxd matrix components [T^]. 

C(f%,Npr) = --±^N[tr{N^ v n vll ) + d(n-2)A b ] 

~ lib g(n ~ 2)ylz + £ ™(g^,X, iM- • ■)» (5) 

where A b ra—A z so that the total yl matches astronomical measurements [23] 

yl = yl 6 + yl z wlO- 56 cm- 2 , (6) 
and the vector potential is defined to be 

A v = r { l a] /[{n-l)ij2A b ]. (7) 
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The "physical" metric g v>l and the fields g ufi , W , g and N are defined by 

gg"" =NN^\ gg u ^ = V~g (I g y ^-h^), tr(h v>t ) = 0, (8) 
g=(±det(g vii )) 1/2d , N=(±det(N Ufl )) 1/2d , + for even d, - for odd d. (9) 

Note that (HJ) defines g M " unambiguously because g = [±det(gg ,l ' / )] 1 / d ( n ~ 2 \ 
The symmetric metric g vp is used for measuring space-time intervals, covariant 
derivatives, and for raising and lowering indices. The C m term is not to include 
a tr(F pa g all g pu Fv f j i ) term but may contain source terms with the usual coupling 
to A v and <7„ M . Tensor indices are assumed to have dimension n=4, but as with 
the matrix dimension "d" , we will retain "n" in the equations to show how 
easily the theory can be generalized. The non-Abclian Ricci tensor in |5]) is 
chosen to have special symmetry properties to be discussed later, 

pr j^P 

'tS pa _ pa _r_ 1_ per riQ _■_ 1 pa per pier pia _ \ rv \ [PfA /-i r>\ 

'^ufi ^ ufi.a ^(0(1/),^)+ 2 (™)+ 2 (™) va 1 ff/j / _ j\ ■ l lu J 

For Abelian fields the third and fourth terms are the same, and this tensor 
reduces to the Abelian version in [T3] . This tensor reduces to the ordinary Ricci 
tensor for , = and r£, , =0, as occurs in ordinary general relativity. 

The determinants g = det(g„ M ) and N = det(N utl ) are defined as usual 
but where N ufl and g are taken to be nd x nd matrices. The inverse of 
N vfl is defined to be N^^ kvl = (l/N)dN/dN^ t k where i,k are matrix indices, 
or N^ LV = (l/N)dN/dN ufl using matrix notation. The field N^ v satisfies 
the relation N^ km N virTj = S£S k , or N^ U N V(J = 5%I using matrix notation. 

Likewise g V(J is the inverse of g M " such that g^g,^ = S£I. Assuming N aT — 
T^N vfJi T^ for some coordinate transformation T£ = dx u / dx a , the transformed 
determinant N = det(N aT ) will contain d times as many T£ factors as it would 
if N aT had no matrix components, so N and g are scalar densities of weight 
2d. The factors N and g from ((9]) are used in (|5|) instead of \/—N and \J—g to 
make the Lagrangian density a scalar density of weight 1 as required. 
For our theory the Maxwell- Yang-Mills field tensor f v ^ is defined by 

g/"M = iNN^A\ /2 /V2. (11) 

Then from {§]), g pv and are parts of a total field, 

(N/g)N^ p = g^ + if^V2Ar b 1/2 . (12) 

We will see that the field equations require fi/ fl ^2A[ IJ , jU ] + i\J2Ab [A v , A^} to 
a very high precision. From ([2|) this agrees with Einstein-Maxwell- Yang-Mills 
theory when \f2A~b — g c /Ti\/2d. Using d=2, g c — e/sin9 w and ^ gives 



where l P = ^/Gh/c 3 = 1.616 x lCT 33 cm, a = e 2 /Tic= 1/137 and sin 2 W = .23. 
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It is helpful to decompose into a new connection J 1 " , and A v from (J7|, 

= + {8«A V - 62A„) iy/2A b , (14) 
where = (Sffi^ - ^f^)/(n-l). (15) 

By contracting (fTSj) on the right and left we see that has the symmetry 

^va—^(la)—^avi (16) 

so it has only n 3 — n independent components whereas r"^ has n 3 . Substituting 
the decomposition (fT4)) into (fTOjl using (jllOp from Appendix lAl 

K Ufl (r) = K Vfl (r) + 2A M i^/2A b + 2A b {A v ,A ll ] 

+ ([A a ,r^] ~ [A {v ,r^ a ])i^2A h . (17) 

Using (fT7|) , the Lagrangian density (O can be rewritten in terms of jT"^ and 
Aa from (fT5T7|) . 



£ = — iV 

167rd 



ir(AT H ^(7^ M + *V2A+ 2A b [A v ,A ll ] 

+ ([A a ,r^] - r° )o ]) *a/2^ )) + d(n-2)A b 



1 

16tt 1 



;(n-2)4, + £ m (g Miy ,X,^,0. • ■)• (18) 



Here 71^^ — TZ utl (r), and from (TITj)) our non-Abelian Ricci tensor (TIT))) re- 
duces to 

•7-7 _ r^a I 1_ j-\a j-*a < _1_ j^a j^a j^a -pa /-i q\ 

From (|14ll6p . and _4„ fully parameterize and can be treated as in- 
dependent variables. The fields AW H (^) and AW H M (or g w " and /" ") fully 
parameterize N vfl and can also be treated as independent variables. It is sim- 
pler to calculate the field equations by setting SL/ST"^ = 0, 5C/8A V = 0, 

6C/5{NN<^) = and SC/S(NN^^) = instead of setting SC/ST^Q and 
8C/SN vfi = 0, so we will follow this method. 



3 Invariance properties of the Lagrangian density 

Here we show that the Lagrangian density is real (invariant under complex con- 
jugation), and is also invariant under U(l) and SU(d) gauge transformations. 
The Abelian Lambda-renormalized Einstein- Schrodinger theory comes in two 
versions, one where F?^ and N vll are real, and one where they are Hermitian. 
The non-Abelian theory also comes in two versions, one where 2"^„ and N VIX 
are real, and one where they have nd x nd Hermitian symmetry, r^* k = r" kvi 



6 



and N* ijlk = N^kvi, where i, k are matrix indices. Using matrix notation these 
symmetries become 

pa*_paT pa*_paT iu« njT TyH^ii; * __ l\H^M T (OQ\ 

where "T" indicates matrix transpose (not transpose over tensor indices). We 
will assume this Hermitian case because it results from A z < 0, At, > as in 
(fT3f . From (|20I8I11I7|) the physical fields are all composed of d x d Hermitian 
matrices, 

V/J.T * T fun* muhT f* _ rT pot* potT a* jT 

Hermitian f v)Ji and A v are just what we need to approximate Einstein-Maxwell- 
Yang-Mills theory. And of course g 1 ^ and g ufl will be Hermitian if we assume 
the special case where they are multiples of the identity matrix. Writing the 
symmetries as N* ifik = , g* t>A = g ukfli = g flkv i, and usin S the result that 
the determinant of a Hermitian matrix is real, we see that the ndx nd matrix 
determinants are real 

N* = N, g*=g, g*=g. (22) 

Also, using ([2U)l and the identity MfM^f = (M2Mi) T we can deduce a re- 
markable property of our non-Abelian Ricci tensor (|10[) . which is that it has 
the same ndxnd Hermitian symmetry as -T"^ and N^, 

(23) 

From the properties (j23l2QI22j) and the identities tr(MiM 2 ) = tr(M 2 Mi), 
tr(M T ) = tr(M) we see that our Lagrangian density or (fT5| is real. 

With an SU(d) gauge transformation we assume a transformation matrix 
U that is special (det(U) = 1) and unitary (WU — I). Taking into account 
(|3l7ll4p , we assume that under an SU(d) gauge transformation the fields trans- 
form as follows, 

(24) 

(25) 

(26) 
(27) 
(28) 
(29) 
(30) 

U UfuufJ' 1 , (31) 
-> Uf^U' 1 . (32) 





-> UTaAlU- 1 - 




A-u ~ 


-» UAuU- 1 + 




A - 






na 


- erf ^tr 1 + 


2^1/^jC/- 1 , 


pQ 

1 M ~ 






pen 

1 [an] - 




f (n-l)^- 1 , 




+ u%u-\ 




AW - 


-> UN^U- 1 , 
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N 32 
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Under a U(l) gauge transformation all of the fields are unchanged except 

A v ^Al+-±=<p„ (33) 

V *> A h 

A v ^> A v + -^=f, v , (34) 

r^^r^-2i/^], (35) 

^l-^]-^- 1 )^- (36) 
Writing the SU(d) gauge transformation (f3"Tj) as 

^ = I [/ ] I N 2 l JVa! ^22 ^23 I I C/- 1 ' (37) 

and using the identity det(MiM2)—det(Mi)det(M 2 ), we see that the ndx nd 
matrix determinants are invariant under an SU(d) gauge transformation, 

N^N, g^g, g^g. (38) 

Another remarkable property of our non-Abelian Ricci tensor (j 10)) is that it 
transforms like N Vfl under an SU(d) gauge transformation l[2"7|). as in (|106|) of 
Appendix [A] 

1l Ufl (Ur ^ T U- 1 + 2S^ p U tT] U- 1 ) = UTl VIJ ,{T^ T )U- 1 for any matrix U(x a ). (39) 

The results (|38I39[) actually apply for any invertible matrix U, and do not 
require that det(U) = l or U'U — I. Using the special case U = Ie~' lv in (|39|) 
we see that our non-Abelian Ricci tensor (fT0| is also invariant under a U(l) 
gauge transformation, 

H Ufl (r«-2iI5% ) <p tT] )=n vfl (f«) for any y^). (40) 

From (|39I31I38I40P and the identity tr(MiM 2 ) = tr(M 2 M 1 ) we see that our 
Lagrangian density (O or (fT8|) is invariant under both [/(l) and SU(d) gauge 
transformations, thus satisfying an important requirement to approximate 
Einstein-Maxwell-Yang-Mills theory. 

One of the motivations for this theory is that the A z = 0, C m = version 
can be derived from a purely afhne Lagrangian density as well as a Palatini 
Lagrangian density, the same as with the Abelian theory[5]. The purely affine 
Lagrangian density is 

C{f % ) = [±det(N utl )] 1/2d , (41) 
where N v ^ is simply defined to be 



N Vfl = -H v JA b . 



(42) 



Considering that N^ u = {l/N)dN/dN Ufi , we see that setting 5£/8r^ = 
gives the same result obtained from (0 with A z = 0, C m — 0, 

tr[N^"STl VIM /Sr^} = 0. (43) 

Since gl]) depends only on f^, there arc no 8C/5(NN^ U ) = field equa- 
tions. However, the definition (g2J) exactly matches the 5C/5(NN^ L ") = field 
equations obtained from with A z — 0, £ m = 0. 

Note that there are other definitions of N and g which would make the 
Lagrangian density (|5|) real and gauge invariant, for example we could have 
defined N = tr(det(N u ^)) or N — Det(det(N l/M )), where det() is done only 
over the tensor indices. However, with these definitions the field N^^ v = 
{l/N)dN/dN vll would not be a matrix inverse such that N^ a,/ N l , tl = 6°I. Cal- 
culations would be very unwieldy in a theory where — (l/N)dN/dNv^ 
appeared in the field equations but was not a genuine inverse of N ufl . In addi- 
tion, it would be impossible to derive the A z = 0, C m = version of the theory 
from a purely affine Lagrangian density, thus removing a motivation for the 
theory. Note that we also cannot use the definition TV = det(tr(N„ M )) as in 
[20] because det(tr(N„ p )) and det(tr(7?. !yAt )) would not depend on the traceless 
part of the fields. 



4 The case h Ufl — with nonsymmetric fields 

Let us consider the theory for the special case = 0, or more precisely for 

f^=tr(f^)I/d, ?» = W)I/d. (44) 

In this case A v and NN^"^ are the only independent variables in (p~8|) which 
are not just multiples of the identity matrix /. This assumption is both co- 
ordinate independent and gauge independent, considering (|30|32p . We assume 
this special case because it gives us Einstein-Maxwell- Yang-Mills theory, and 
because it greatly simplifies the theory. With the assumption ([44")) we also have 
= tr{H v ^)I/d, and the term {[A a ,f^] - [A {[J1 f« a }) iy/tt~ b vanishes in 
the Lagrangian density (I18p . It is important to emphasize that any solution of 
the restricted theory ((44)) will also be a solution of the more general theory. 

Setting 8C/SA T = and using the definition (fTTj) of gives the ordinary 
Maxwell- Yang-Mills equivalent of Ampere's law, 

(gf " iV^hiir^A^} = 47r£f , (45) 



where the source current j T is defined by 
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Setting S£/5r^ p — Q using a Lagrange multiplier term tr[fl p r^ p ^\ to enforce 
the symmetry (|T6|) . and using the result tr[(gf UJT ) i u ] = 4ngtr[j T ] derived from 
(|45|3|4p gives the connection equations, 

tr[{NN^ pT \ p + r^pNN^f" 7 + r^ a NN~ lcrT - r$ a NN ApT ] 

^%itr\j^. (47) 



{n-l)A\l 



Setting 6C/S(NN^^) = using the identities N = [±det{NN^ u )} l / d{ - n -^ 
and g= [ide^AW -1 ^"^)] 1 /^" -2 ) gives our equivalent of the Einstein equations, 

2 tr [K(vn) + A b N M + A z g vp ] = &irtr[S vli ], (48) 
where S v ^ is defined by 



( 49 ) 



Setting S£/S(NN^^) = using the identities N = [±det{NN^)] 1 /< n -^ 
and g=[±det(NN<^)] 1 / d( - n - 2 ^ gives, 

+ 5U[„ fM] i v / 2^+ 2yl 6 [A,^] + ^fcA^j = 0. (50) 

Note that the antisymmetric field equations (J3UJ) lack a source term because 
C m in CEE]) contains only gg^ = jVAT 1 ^ f rom flgj^ an d not A/ATM. The 
trace operations in (I47I48|) occur because we are assuming the special case 
(|44| . The off-diagonal matrix components of 5£/8fP p and 5C/5(NN^ P ^) van- 
ish because with (|44|1 , the Lagrangian density contains no off-diagonal matrix 
components of r@ and NN^^"'. The trace operation sums up the contri- 
butions from the diagonal matrix components of /f p and NN^^^ because 
(HI]) means that for a given set of tensor indices, all of the diagonal matrix 
components are really the same variable. 

To put (|45H50p into a form which looks more like the ordinary Einstein- 
Maxwell- Yang-Mills field equations we need to do some preliminary calcula- 
tions. The definitions (|8I11[) of g up and f vp can be inverted to give N up in terms 
of g vp and f vii . An expansion in powers of A^ 1 is derived in Appendix [Bl 

N M =g„ - 2 ~ 2(^)^ trU T P) ) A * 1 + {f)AbV2 • ■ ' (51) 

N Wp] =f„ p iV2A- b 1/2 + (f)A- 1 .... (52) 

Here {P)A^ /2 and (/ 2 ) V are terms like faf^f^A^ 2 and f [v f Aa /^\ 
Because of the assumption (j44[) and the trace operation in (j47|) . the con- 
nection equations (j47|) are the same as with the Abelian theory [13] but with 
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the substitution of tr[f utl ]/d and tr[j v ]/d instead of f vjl and j". Therefore the 
solution of the connection equations from |13j can again be abbreviated as 

A"m) = + V'f)^ 1 ■ ■ ■ f M = (/W ■ • • > (53) 
where is the Christoffel connection, 

r?K = \ 9 aa (g^, v + g^p - g^a). (54) 

Substituting (|53p using (fT02|) shows that as in 13] . the non-Abelian Ricci 
tensor (JTUJ) can again be abbreviated as 

= ^+ (/7')V + • • • , *M = (/")^ 1/2 • • • , (55) 

where = (J 1 ) is the ordinary Ricci tensor. Here (f'f)A^ 1 , [ff")A^ 

and (/")^ 1/2 refer to terms like tr(f%, a )tr(f a Ka )A^\ tr(/°r)M/W);«0V 

and MW]; tt K^- 

Combining (|51I55I6[) with the symmetric field equations (|4"B"|) and their 
contraction gives 

G vll = 8irtr(T ull ) + 2 ^ — ' -3^ — 



d 4""" d 

+A g - l) 5 „ M + (/ 3 )^ 1/2 + (/Y)V + (//") V ■ ■ • , (56) 

where the Einstein tensor and energy-momentum tensor are 

Gj/ju = -Riy^i — ~^9vijl R a i — S UfJ ,— — g u ^S a . (57) 

Here (f 3 )A~ 1/2 , (/'/') V and <Jf")A? ara terms like tr(/'„ f^U),,)^ 2 , 
frC/^aM/'VM* 1 and tr(f ar )tr{f T{ „. „j ;a )^ 1 . This shows that the Ein- 
stein equations (|56[) match those of Einstein-Maxwell-Yang-Mills theory except 
for extra terms which will be very small relative to the leading order terms 
because of the large value Af, ~ 10 63 cm -2 from (fT3"]) . 

Combining (|52I55|) with the antisymmetric field equations ([50)1 gives 

U = 2A M + iV2Jk[A u ,A ll ] + (f)A- 1/2 + .... (58) 

Here (f 2 )A- 1/2 and are terms like f^J^ 1 and MW];") 

From (fT3|) we see that the in Ampere's law (j45|) matches the Maxwell- 
Yang-Mills tensor ((2]) except for extra terms which will be very small relative 
to the leading order terms because of the large value 10 63 cm~ 2 from (I13| . 

Let us do a quantitative comparison of the h vl1 — case to Einstein- 
Maxwell- Yang-Mills theory. To do this we will consider the magnitude of 
the extra terms in the Einstein equations and the Maxwell- Yang-Mills field 
tensor for worst-case field strengths and rates-of-change accessible to mea- 
surement, and compare these to the ordinary terms. In particular we will 
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evaluate extra terms in the Einstein equations ([56")) like tr{f p rj f l7 ( p f v ) p )A b ' , 
M/° u-,a)tr(f a fj,;<r)A^ and tr{f aT )tr{f T ^ u . tl y ia )A^ 1 and compare these to the 
ordinary Maxwell- Yang-Mills term. Likewise we will evaluate extra terms in 
the Maxwell- Yang-Mills field tensor (J58j like f^f^A^ 2 and tr{f [v ^ a ].^)A^ 1 
and compare these to which appears in Ampere's law (|45[) . 

We assume that the worst-case field strengths and rates of change accessible 
to measurement will be purely electromagnetic fields. Also, because we will just 
be doing order-of-magnitude calculations, we will neglect mixing in /^ t „ and we 
will use the electromagnetic coupling constant. In geometrized units with the 
Heaviside-Lorentz convention an elementary charge has e = 1.38 x 10 _34 cto. If 
we assume that charged particles retain f 1 o^e/r 2 down to the smallest radii 
probed by particle physics experiments (10~ 17 cm) we have from (|13l) . 

l/^j^ 1 / 2 ~ ^ 1/2 e/(lQ- 17 ) 2 ~ 10- 31 , (59) 
l/Vi/Al 2 ^ 1 ~ 4^V(10- 17 ) 2 ~ 1CT 29 , (60) 
If^i/fW^ 1 ~ 6AjV(10- 17 ) a ~ 1CT 29 . (61) 

The fields at 10~ 17 cm from an elementary charge would be larger than near 
any macroscopic charged object. Here / o is assumed to be in some standard 
spherical or cartesian coordinate system. If an equation has a tensor term 
which can be neglected in one coordinate system, it can be neglected in any 
coordinate system, so it is only necessary to prove it in one coordinate system. 
So for electric monopole fields, the extra terms in the Einstein equations (|56|) 
must be < 10 -29 of the ordinary Maxwell- Yang-Mills term. Similarly the extra 
terms in the Maxwell- Yang-Mills field tensor ([55)) must be < 10 -29 of f vfJ ,. 
Also, for the highest energy electromagnetic waves known in nature (10 20 eV, 
10 34 Hz) we have from Iflgjl. 

I/V/AI'V ~ {E/hcfA- 1 ~ 1(T 13 , (62) 
l/W/Z'olV ~ (E/hcfA, 1 ~ 10- 13 . (63) 

So for electromagnetic waves, the extra terms in the Einstein equations (|56p 
must be < 10 -13 of the ordinary Maxwell- Yang-Mills term. Similarly the extra 
terms in the Maxwell- Yang-Mills field tensor ([58)) must be < 10 -13 of 
which appears in Ampere's law f|45[) . 

From this analysis we see that these extra terms in the field equations 
(I56I58I45[) are far below the level that could be detected by experiment for 
worst-case field strengths and rates of change accessible to measurement. At 
least we have made great efforts to find an experiment in which these extra 
terms would be evident, and we have been unable to find such an experiment. 
As shown in [13] . the ordinary Lorentz force equation can be derived from 
the divergence of the Einstein equations for the purely electromagnetic case 
of this theory. In [13) we also presented an exact electromagnetic plane-wave 
solution which is identical to its counterpart in Einstein-Maxwell theory. And 
in [13j we presented an exact electric monopole solution which matches the 
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Reissner-Nordstrom solution except for additional terms which are < 10~ 65 of 
the usual terms for worst-case radii accessible to measurement. 

We wish to emphasize that the C m term in ([5]) allows coupling to addi- 
tional fields via a symmetric metric <? M „ and Hermitian vector potential A^, 
just as in Einstein-Maxwell- Yang-Mills theory. Our £ m can contain the same 
fermion and Higgs field terms as in Weinberg-Salam theory or flipped SU(5) 
GUT theory. And when we do this we will get the same energy-momentum 
tensor (|49|57|) in the Einstein equations (l56|) . and the same current (|46|) in 
the Maxwell- Yang-Mills equivalent of Ampere's law (|45I58|) . In addition, the 
equations of motion of fermion and Higgs fields will be unchanged, and the 
components of A^ will mix and acquire mass in the usual way (the A^ mass 
terms will get lumped into j T in (|45|46| ). 

One aspect of this theory which might differ from Einstein-Maxwell- Yang- 
Mills theory is discussed in detail at the end of section 5 of [13] for the purely 
electromagnetic case, although it is unclear whether it is really a difference or 
not. To see what this is we take the curl of ([58]) . in which case the 2Ay fl ^ v \ term 
falls out, and from the $ vyL and (f")A^ terms we get[13j. 

f[vn, a ] = {-f[un, a ); a ';c + apparently negligible terms)/2yl 6 .... (64) 

This is similar to the Proca equation with the field 6 T = e TV ^ a /r 1 , J[iia ]/4. It 
suggests that the theory may allow 9 T Proca waves with mass from (|64I13|) 
close to the Planck mass. For d = 2 and g c = e/sin6 w we get 

1 9c fa „ /r t __ 1 1 w1 . 18/ 



uproca = v 2A b — \hr-,, M Proca = hujp rooa = l.l XlO GeV. (65) 

lp e V 2d 

Using a Newman- Penrose 1 jr expansion of the field equations we have shown 
that continuous-wave solutions like 9 T w e T sin{kr — ut) / r do not exist in the 
theory [19], but it is still possible that wave-packet solutions could exist. If 
wave-packet 9 T solutions do occur, a calculation in [T3] also suggests that they 
might have negative energy, although this calculation is really based on the 
assumption that 9 T ~ e T sin(kr — uit) / V solutions exist, and some questionable 
assumptions about terms being negligible. If wave-packet 8 T solutions do exist, 
and if they do have negative energy, there is still a possible interpretation of 
the 9 T field as a built-in Pauli-Villars field, with a cutoff mass (|65|) which is 
close to Mpi anc k = 1.22 x 10 19 GeV commonly assumed for this purpose. 

The additional cosmological constant A z in our Lagrangian density (|5I18[) 
could have several contributions. If there was a contribution from zero-point 
fluctuations it would be approximately 14. 15. 16,17] 

_ ^V-p ( fermion _ boson \ /gg% 
z0 ~ 2tt \ s pi n states spin states / ' ^ ' 

where lj c is a cutoff frequency and I p = (Planck length) . Assuming the Pauli- 
Villars ghost idea discussed above, Lo c — uop roca from (|6"5]) , A z fa —At from ([2 
d = 2, g c = e/ sin9 Wl and the particles of the Standard Model gives 

A z0 _ Up roc Jp , nc OQ , _ a fg c \ 2 



A z 7r 2dn V e 



(96- 28) = — — (96 -28) = .17. (67) 
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So by this calculation, zero-point fluctuations would only contribute about 
17% of A z . Additional contributions to A z could perhaps come from Higgs 
field vacuum energy and additional unknown fields. It is unclear how this 
calculation would work out for flipped SU(5) GUT theory. Note that ujp roca 
and At, from (f65|) depend on the coupling constant g c , so their values should 
"run" with frequency. The contribution to A z from zero-point fluctuations (|67[) 
would be slightly modified if we used a "bare" g c calculated at top roca instead 
of a low energy value. Also note that the Pauli-Villars ghost idea might not 
be necessary or correct, in which case we could make (p7|) closer to 100% by 
assuming a slightly larger u> c . 



5 The case \h vtl \ <§C 1 with symmetric fields 



This theory definitely differs from Einstein-Maxwell- Yang-Mills theory in that 
the symmetric fields can be non-Abelian, with traceless components. To in- 
vestigate this let us calculate the field equations with the Lagrangian den- 
sity (jT8TT9l> and the special case A v = 0, = 0, = 0. Setting 
SC/5(NN^^) = and using N = g = [±det{NN^^)] 1 / d ^ 1 ^ gives our 
equivalent of the Einstein equations, 

1 



where S vp 



5 C m 



„ SCr, 



(68) 
(69) 



For present purposes we assume S vfJi — and yl = 0. Setting SL/ST^ =0 using 
,=0, N^v] = Q 7 A u =0 and © gives the connection equations [15]. 



t, 



TV Tip 
I 1 V0 

rr$ a = o. 



(70) 



We will only consider the case where the traceless components are small, 
similar to linearized gravity, 



|/^|«i, |fl?J<|r»|, f"=ir"+H"+o(h 2 ). 



(71) 



Here h vp is defined in (|8|) and is the Christoffel connection (f54|) formed 
from the physical metric g vpj with no traceless components. The connection 
equations ([70)1 to C(/i) are 



1 pT 

'2 ^ 
1 



-r: 



/3 a 1 



2 ^ 

1 TTT 

2~ 



/l.T_ 



99 



hp _ 



2 v y up 2 



"0 



/3a 1 



\j—gh pT r$ a + I^fl^fl^ - 0. (72) 
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Using (V-fl),/3 = \F~9 r p a and dividing by ^g gives 

= hT. p - - gTKfi + 9 pT H$ a . (73) 

Combining the permutations of this gives 

= (/i W A:/3 — H^px — H\up + g^\Hp a ) 

— {hp^-x — Hpxu — H^px + gpijHxa) 

— (h\/3;w — Hxuip — Hpxuj + gxpH^a) (74) 
= + h^x-p — hp^.x — hx/3; U + 9u>\Hp a — gp^Hxa ~ 9\pH" a . (75) 

Contracting this with g^ gives 

= 2HX u -K. x -nH% a => fl^ = _l_j#. A . (76) 

So the 0(/i) solution of the connection equations ((70)) is 

H a v P — ^p l av^^^^ia-,v~~^ l v^\a) + i.9av^>u;n ~^~9p.a ^lo^u ~ 9t>fJ^>u;a) ' (77) 

Assuming Z 1 " = J T"+H"+K" +0(h 3 ) and using a similar method [TH] gives 
the 0(h 2 ) solution of the connection equations (|70p. 



^/3rp = ^ [-9pt(KK),P + 90p{KK),T + 9rp{KK),P 

+ {hpa-T— h aT .p)hp + hp(hp a . T — h aT .fj) + h pT .^h a p + hph pT . a ]. (78) 

The field equations for h vp are found by substituting the 0(h) solution (|77|) 
into the traceless part of the exact field equations (168)) and using (|102I76I8|) 



= 2 [H^. a - H« M ] (79) 

s a 4- 

(n-2)" 



-h UIL -a; a + 2/l a ( I/;jli ) ; Q + . ^ oA ffi/ M ^-r;a; a - (80) 



Contracting this equation gives 

K. a . a = (2-n)hl ;T . a . (81) 
So we can also write the field equations as 

= —hv^a™ + 2ha(v;n); a — 9vph T a - T a ■ (82) 

Now let us assume that we can ignore the difference between covariant 
derivative and ordinary derivative. In that case (|80I82|) match the "gauge in- 
dependent" field equations [55] of linearized gravity, but with a non-Abelian 
h vp . In linearized gravity one often assumes the Lorentz gauge 

K a; a =0. (83) 
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Here we do not have the same freedom because in the coordinate transforma- 
tion x v — > x u ', /i^— > h Ufl — £v tP ,— £n;v the parameter cannot be traceless 
like hvfj,- However, we can still seek solutions which satisfy (f83|) . Analogous 
with linearized gravity we have a z-directed plane-wave solution 



l l> P 



sin(ujt — kz) 







h A 
h> 

Vo 







k = a;, 



(84) 



and a static spherically symmetric solution 



/4M/r 0* 





V 0, 



(85) 



Here h + ,h x ,M are constant traceless Hermitian matrices. Note that (1551 vi- 
olates Ift^l^Cf near r — 0, so a corresponding exact solution may not exist. 

To find an effective energy-momentum tensor for h vp we extract the 0(h 2 ) 
componentspjj] from the exact field equations flBSJ) using (|f 02I77I76I78I8|) . 

8nS Tp = -tr[K%. a - K a a(r , p) + K P K<* - K<*K P \ (86) 
= tr[-g^(Kh»). ia . t a /4(n-2) + g pT {h^ a K). a /2{n~2) 

+ h Tp °h a a ,j2{n-2) - hZ. p K. T /4(n-2) + h^ a h%j2 - K. a K p »/2 
+ K ;T h« p /4}. (87) 

So the effective energy- momentum tensor is|19] 

8nf Tp = 8^ (S TP - \gr P S^ (88) 
= tr[-g pr (KK). a ., a /8 + g pT (Kh^), a , a /8{n-2) + g pr (h^K), a /2 
- gpT h^h%j8{n-2) - g pT h^, a h%j4 + g pT h^. a h^. a /8 

+ h rp ?h%j2{n-2) - h u u . p K. T /4(n-2) + K ;a h%j2 - K. a K m a /2 
+ h°. T h% p /4\. 

From the field equations ([50]) we get 



= *r[(-^;% + 2^. p) . a + g pu hZ ;a . a /(n-2))K] (90) 

- ir[-^ ; a ^ + K. p h v T + /i«X + hZ.y pT /{n-2)], a 

- tr[-h up . a K. a + K, p % a + h%»K; a + K, a Kr;J(.n-2)}. (91) 
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Using tr(M\M2) = tr{M%M\), the symmetrization and contraction of (|91l) arc 
= tr[-(h vp K), a /2 + K. {p h" T) + hffi );v + h». a h pr /(n-2)] ;a 

- tr[-K p . a K;a + K.{ P K); a + K>r; a + K,° V ; a/(»" 2 )1 > ( 92 ) 

= tr["( W);72 + 2^ ; >- + (/&/£); a /2(n-2)] ;a 

- tr[-K. a K ;a + 2h a v .^°, a + ^ ; <X ;a /(n-2)]. (93) 

Adding to ([89]) the expression (j92j) /2— a^ (f9"3"]) /8 gives a simpler form of the 
effective energy- momentum tensor [T9] which is valid when S^^— in (I68|) . 

8^f rp = ir[-. 9pr (^^) ;a; 7l6 + g pT (h^h^), a /4 + g pT (h^h a a ). a . a /16(n-2) 

~ (KafrK))" + («);a;74 " ( WO"/ 2 + (h^h^);^ 

+ KXJ^ ~ KpKrKri-2) + K M , a h\ T) /2\. (94) 

Averaging over space or time, covariant derivatives commute and gradients do 
not contribute 26 , so the averaged effective energy- momentum tensor is 

^<f Tp >=<tr[K,^jA-K^ (95) 

This result is the same as for gravitational waves [2"B] but with a non-Abelian 
h vp . From (|95I83I4|) we see that the solution ([84]) has positive energy density, 

8tt < f 00 > - < tr[hZ. h° ] > /4 (96) 

= < tr[h\. hl. + hl. hl. Q + h\,Ji\, + /i2; ^ ; o]/4 > (97) 

= tr\h\ + h 2 x ]uj 2 /A > 0. (98) 

While solutions like (|84I85|) have not been observed, one must remember 
that gravitational waves and black holes have not been observed directly either. 
Solutions like (|84l85p do not rule out the theory. In fact if there is an exact 
solution corresponding to (I85p . it might be a possible dark matter candidate. 

6 Conclusions 

The Einstein-Schrodinger theory is modified to include a cosmological con- 
stant A z which multiplies the symmetric metric, and by allowing the fields to 
be composed of Hermitian matrices. The additional cosmological constant is 
assumed to be nearly cancelled by Schrodinger's "bare" cosmological constant 
Ab which multiplies the nonsymmetric fundamental tensor, such that the total 
"physical" cosmological constant A = Ab + A z matches measurement. If the 
symmetric part of the fields is assumed to be a multiple of the identity matrix, 
the theory closely approximates Einstein-Maxwell- Yang-Mills theory. The ex- 
tra terms in the field equations all contain the large constant Aj, ~ 10 63 cm -2 
in the denominator, and as a result these terms are < 10 -13 of the usual 
terms for worst-case fields and rates of change accessible to measurement. Like 
Einstein-Maxwell- Yang-Mills theory, our theory is invariant under 17(1) and 
SU(d) gauge transformations, and can be coupled to additional fields using a 
symmetric metric and Hermitian vector potential. 
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A Some properties of the non-Abelian Ricci tensor 

Substituting = + 7^ into JT9J) gives 

77 ( p\ — pa _ pa , 1_ per pa , 1 pa pa _ ptr pa fqq\ 
l\-v[± K 1 } — 1 v\i,a 1 a(v,{j.) ^ 2 ^ 2 vot (T fl V yy I 
I pa _i_ 'yet \ / pa i y*a \ / pa i ■v^a *i / pa i 'ya \ 

■|(^ + ^ M )(^+T^) + I 



+ r^)(r? a + K a ) + kr? a + r»j(r^ + j~) (100) 



2 



l vfA,a a(v,fi) vat J er/j, ± i/a. i a}i va cr [i. 

V v/i era. ' * vfj,* a a. ' l>h <ra ' era u ~t~ ^ aa i/fi \ aa.*- v{i) \^^) 



7~> / 7-i\ i yQ "y^Q "^(7 "fGL I J- -yCT ■ J_ -yd "yfT /-I f>Q\ 

— -Tti/^^ ; -T J-vf_ L - a — 1 a (u;^) ~~ 1 va 1 <7fi "T ~ 1 i/fi I aoi > ~ era 1 Vfj,- \ 1[JZ ) 

Substituting the SU(d) gauge transformation -» = l/f^JJ- 1 + 2£«£/ jM ] C/" 1 

from J27D into 7^.^ proves the result H39D , and the result J40I I for a U(l) gauge transformation 
Fp T -* r* T -2U5™ p v jT] follows for the special case U = Te~*V. 

t pr ipP 

f-y / 1 p\ i pia i pa , 1 t per i no , 1 t p~a ( pa t per 4 pa YjjA [pp] / 1 HQ \ 

Kvpy J- ) — J- Ullia 1 { a {v), t i) 2 (fa)" 1 " 2 Oa) 1 1 va 1 tr^ (n— 1) 

+ l^aa)^" 1 {urz^u- 1 + s^u^u- 1 - s^u- 1 ) 

"(A) ( C/? fr,] C/ " 1 +( rt - 1 ) C/ .- C/ " 1 )( C/? [p fl ] C/_1+(n " 1)C/ -" C/ ~ 1 ) (104) 

rr I pa pa 1 1 pa pa ■ 1 pa pa pa pa L T[/ J LPMJ | rr — 1 

— u I 1 v^,a J (a(i/), M ) "+" 2 "M (o-a) " t " 2 " Q CT M (n-1) J 

- ±u f, a M' 1 - ±uT< a m~} - -u v r, a M' 1 - ±ur, a Mf 

2 2 ( a(y ) 2 2 ( q m) ' 

+ ~U nf, 01 M' 1 - -U v f, a M' 1 

2 ,M ("«) 2 '' ia ^ 

_ Xjipa TT-1 1 krrpa TT-1 
2 U1 (va) U ,p. + 2 W 

+ urZaUj- ur^u- 1 - u, a rz tl u~ 1 + u^r^u^+p-n^uj-u^u- 1 

+ Ur\ Tv] J7- 1 - C/,„f f w] C/- 1 + (n-X)U, v U-* (105) 
= UKvuffiu- 1 . (106) 



18 



Substituting = + (S£A U - <*EMm) from gU into 7e„ M and using f ° 



r^ a) =r^ from JTBJ with the notation [A, B] = AB-BA gives JTTt , 



-r,?„B° 



,.- H [pm] 



( 



+ (<^A,- v 7 ^),- " 

+2A b (n-l)AuA fl 



(107) 



(108) 



r, 



i 



1 



+2^[„, M ]iV 2/1 i> 

-ri a {&lAa-^A„)i^iA h 

~{6° a Av~5ZA a )r«^^/2A b 
+2A b (n-l)AuA fl + 2A b ((2-n)A v A IA - A^Au) 
K„„(f) + 2A[ v , lt ] i<s/2Jk + 2Ab[Au, Ap] 

+([A»,f? / J - lA (l ,,r^ a ])iy/2A b . 



(109) 
(110) 



B Approximate solution for N vfJ , in terms of g„ M and /„ M 

Here we invert the definitions (181111 ) of g y)J and fvji to obtain II51I52I I. the approximation of 
Nvfi in terms of and f vfl . First let us define the notation 

f"v = fH\/2A' b ~ 1/2 . (ill) 

We assume that \ f n\ 1 for all components of the unitless field f" and find a solution 
in the form of a power series expansion in f" u,- 

For the following calculations we will treat the fields as nd X nd matrices but we will 
only show the tensor indices explicitly. Lowering an index on the right side of the equation 
(N /g)N^^=g^ + f^ from (T3) we get 

(N/g)N^ a =5%I-f> a . (112) 

Using f a a =Q, the well known formula det(e M ) = exp (tr(M)), and the power series ln(l — 
x) = -x - x 2 /2 - x 3 /3 ... we get [27). 

ln(det(I-f)) = tr(ln(I-f)) = -Ur{f"J%) + (f 3 )... (113) 
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Here the notation (/ 3 ) refers to terms like tr(f T a f a a f t)- Taking ln(det()) on both sides 
of d!12t using the result (1113 1 . the definitions |(9j, and the identities det(sM)= s nd det(M) 
and det(M~ 1 ) = l/det(M) gives 

ln(det[(N/g)N^ a ]) = ln((N / 'g)"^ 1 ) = -\ tr(/> CT /%) + (/ 3 ) . . . , (114) 

I«[WS)1 = , tr(f P «f%) + if 3 ) ■■■■ (H5) 

2d(n — 2) 

Taking e 31 on both sides of this and using = 1 + x + x 2 /2 . . . gives 

= ^ * r (/ p CT /%) + (/ 3 ) • • ■ ■ (us) 

2d(n — 2) 

Using the power series (1 — a;) - 1 = 1 + x + x 2 + x 3 . . ., or multiplying by l|112|l on the right 
we can calculate the inverse of II112H to get [27] 

(g/N)N% = 5" I + /% + f v J\ + (/ 3 ) ■ ■ • • (H7) 
Lowering this on the left gives, 

*fup = (A7g)(g^ + fvy. + jW*% + (/ 3 ) ■ ■ •)■ (H8) 
Here (/ 3 ) refers to terms like f va f a CT f" M . Using (1441 1 181 1161 1 1 1 I t we get the result 1511521 . 
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